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1 Introduction 

QCD at non-zero temperature and density is a rich subject discussed in many 
talks at this meeting: from the introductory lesson describing the theoretical 
discovery of the quark gluon plasma El to the^ subtle interrelations of gauge 
fields dynamics and charmonium suppresions □, from the detailed analysis of 
the exotic phases at large baryon density to the fascinating aspects of the 
QCD phase diagram with many flavor from the time-honored theoretical 
subjects of the interplay between topology and confinement □ to the challeng- 
ing aspects of the interrelations of chiral and axial symmetries El. These notes 
mostly aimed at the non-specialist might hopefully provide, besides some in- 
troductory material, a path among general idea as well as among the many 
new results presented in the lively topical session on the QCD phases. 

The material is organised as follows: Section 2 reviews the imaginary 
time formulation of field theory thermodynamics, and gives the functional 
integral representation of the partition function Z. The representation of Z 
naturally leads to a theoretical suggestion: the concept of universality which 
is singled out in a small Section on its own, because it is so important. Sec- 
tion 4 discusses the QCD chiral transition. Section 5 the QCD deconfinement 
transition, and Section 6 how to put together the two, eventually describ- 
ing the complex phenomenology of the real world high temperature QCD 
phase transition. The equation of state is briefly mentioned here. Section 6 
introduces the method used to obtain all of the results of the previous Sec- 
tions: the lattice regularization. We shall see that introducing temperature 
is straightforward while a non-zero density poses specific problems. We will 
only mention the main proposals to circumvent such problems, referring to 
the most recent reviews Ej for details and recent results. 



2 Formulation 



Here we will concern ourselves with the path integral representation of the 

partition function. The basic property of cquilibriiim field theory is that one 
single function Z ( the grand canonical partition function): 

z = z(v,r,Ai) (1) 

determines completely the thermodynamic state of a system according to: 

„ dlnZ , , 

--'-^ (^) 

(4) 

E = -PV + TS + nN (5) 

while physical observables < O > can be computed as 

<0>= TrOp/Z (6) 

Any of the excellent books on statistical field theory and thermodynam- 
ics can provide a more detailed discussion of these points. I would like to 
underscore, very shortly,that the problem is to learn how to represent Z at 
non-zero temperature and baryon density, and how to design a calculational 
scheme. 

2.1 Z at non-zero temperature and chemical potential 

Z is the trace of the density matrix of the system p 

Z = Trp (7) 

p = e(-»-''*)/T (8) 

H is the Hamiltonian, T is the temperature and TV' is any conserved number 
operator. 



2.2 Chemical potential, relativistic and non-relativistic 

It is worth stressing the main differences between the non relativistic and 
relativistic meaning of a chemical potential /x, eqs (5) and (8). 

In a non relativistic setting, the chemical potential tells us the energy 
'cost' of adding an extra particle to the system: we thus have a different 
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chemical potential for each type of particle. The term which is added to the 
Hamiltonian is fJ-kNk, Nk being indeed the number of particles belonging to 
the fc-th species. 

In a relativistic setting, particles can be created and destroyed, thus los- 
ing their individuality, so to speak: it only makes sense to have a chemical 
potential coupled to the 0-th component of a conserved current. 

In practice, we will mostly concern ourselves with the current = ipj^ip, 
i.e. Jo = V'7oV' = V'^V' "will be the density of fermion (baryon) number, 
or, more precisely, the difference between the number of fermions and an- 
tifermions: J Jq — N ~ N . The chemical potential will then appear in the 
Lagrangian as a term /i Jq and we recognize immediately that a change of sign 
for the chemical potential corresponds to an exchange particles-antiparticles. 

A second difference between non-relativistic and relativistic approach con- 
cerns the magnitude itself of /i, which in a relativistic theory contains the 
particle rest mass. 

From a technical point of view, a finite density is easily handled in a non- 
relativistic setting. From a relativistic field theoretic perspectives there are 
instead important technical problems. We will discuss these problems later 
on, and we just anticipate here that the explicit breaking of the particle- 
antiparticle symmetry induced by a chemical potential induces the loss of 
positivity of the effective gluonic Action-i.e. the breakdown of the condition 
which makes lattice calculations possible. 



2.3 Temperature 

Consider the transition amplitude for returning to the original state (j)a after 
a time t: 



{x,t) — (pa (x) 



d(l)e ■fo'^' / d'^{^is,t)S^-H(^,4>)) 

{X.0)=4>a{x) 

(9) 

Compare now the above with expression (2) for Z, and make the trace explicit: 
Z = Tre-Pi"-'^^) = jd^,< </.,|e-''(^-''^)|0, > (10) 
We are naturally lead to the identification 

P^^^it (11) 

We note - anticipating the discussions of Section below - that studying 
nonzero temperature on a lattice c2l is straightforward: one just takes ad- 
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vantage of the finite temporal extent of the lattice, while keeping the space 
directions much larger than any physical scale in the system. 

By introducing the integral S{(j), ■ip) of the Lagrangian density (from now 
on we will always use 1/T as the upper extreme for the time integration) 

5(0, V)= ! ' dt I d^xCi(j),ip) (12) 



Z is written as 

' d(l)d^e-^'^'l'''>'^ (13) 



The only missing ingredient are the boundary conditions for the fields: ba- 
sically, they follow from the (anti)commuting properties of the (fermi)bose 
fields which imply 

^(f,0) = ^(f,/3) (14) 

for the bosons and 

Vi(f,0) = -Vi(f,/3) (15) 

for the fermions. 

Fermions and bosons obey antiperiodic and periodic boundary conditions, 
respectively, in the time direction. 

The expression above, together with the boundary conditions just intro- 
duced, is the key to field theory thermodynamics. 



3 Universality 

It is intuitive that when the smallest significant length scale of the system I » 
1/T the system becomes effectively d-dimensional. Moreover, the description 
of the system can be effectively 'coarse grained', ignoring anything which 
happens on a scale smaller than I . 

This can become true whcin the; system is approaching a continuous tran- 
sition: the correlation length of the system ^ is diverging. In such situation 
all the physics is dominated by long wavelength modes. Not only the system 
gets effectively reduced, but the coarse graining procedure become doable. As 
an effect of this procedure, systems which are very different one from another 
might well be described by the same model, provided that the long range 
physics is regulated by the same global symmetries: this is the idea of uni- 
versality which provides the theoretical framework for the study of the QCD 
transition in two interesting (albeit non physical) cases which we review in 
the next two sections. 



4 



4 QCD chiral treinsition (m^ = 0) 

Let us recall the symmetries of the QCD action with Nf flavors of massless 
quarks, coupled to a SU{Nc) color group: 

SU{N^)c X SU{Nf) X SU{Nf) x Z^iNf) (16) 

SU{Nc) is the gauge color symmetry. SU{Nf) x SU{Nf) x Za^Nj) is the 
flavor chiral symmetry, after the breaking of the classical C/a(1) symmetry to 
the discrete Z^(A^/). 

We want to study the realization and pattern(s) of breaking of the chiral 
symmetries and we would like to know the interrelation of the above with the 
possibility of quark liberation predicted at high temperature and density. 

4..I Ordinary conditions: zero temperature and density 

In normal conditions (zero temperature and density) the SU{Nf)]^ x SU{Nf)]i 
chiral symmetry is spontaneously broken to the diagonal SU{Nf)L+R-ljet us 
note the isomorphism: 

SU{2) X SU{2) = 0(4) (17) 

which shows that the symmetry is the same as the one of an 0(4) ferromagnet. 
The relevant degrees of freedom are the three pions, and the sigma particle, 
and the effective potential is a function of cr^ + |7rp in the chiral space. Once 
a direction in the chiral sphere is selected (say in the a direction) chiral sym- 
metry is spontaneously broken in that direction, according to the equivalent 
patterns: 

SU{2)r X SU{2)l ^ SUi2)L+R (18) 
0(4) ^ 0(3) (19) 

Massless Goldstone particles (in this case, the three pions) appear in the 
direction orthogonal to the one selected by the spontaneous breaking. 

4-2 Increasing T 

Disorder increases with temperature. Then, one picture of the high T QCD 
transition can be drawn by using a ferromagnetic analogy of the chiral tran- 
sition: -ipip can be thought of as a spin field taking values in real space, but 
whose orientation is in the chiral sphere. Chiral symmetry breaking occurs 
when < ■tpip >^ 0, i.e. it corresponds to the ordered phase. By increasing T, 
< ■tptp >— > 0, and the 0(4) symmetry should be restored. 
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Combining this symmetry analysis with the general idea of dimensional 
reduction, Pisarski and Wilczek til proposed that the high temperature tran- 
sition in two flavor QCD should be in the universality class of the 0(4) sigma 
model in three dimensions. At high temperature when symmetry is restored 
there will be just one global minimum for zero value of the fields, and pion 
and sigma become eventually degenerate. 

We have however to keep in mind possible sources of violation of this 
appealing scenario and, all in all, one has to resort to numerical simulations 
to measure the critical exponents, and verify or disprove the 0(4) universality. 
In turn, this gives information on the issues raised for instance in as well as on 
the possible restoration of the axial anomaly, see the discussion in Meggiolaro's 
review. 

In practice, one measures the chiral condensate as a function of the cou- 
pling parameter /3, which in turns determines the temperature of the system. 
This gives the exponent Pmag according to 

< >= B{I3 - pcf"^-^ (20) 

The exponent 5 is extracted from the response at criticality: 

<4)il) >= Am^l^;l3 ^ Pc (21) 

The results for the critical exponents compare favorably with the 0(4) results 
fimag — -38(1), 5 = 4.8(2), and definitively rule out mean fields exponents 
(which would have characterized a weak first order transition). However, the 
results can still be compatible with 0(2) exponents, which would signal the 
persistence of some lattice artifact, and of course it is still possible that the 
final answer do not fit any of the above predictions, for instance if we were 
just observing some crossover phenomenon. 

In conclusion the symmetry analysis of the (two flavor) QCD transition 
gives a definite predictions for the value of the critical exponents, which is 
possibly (slightly) violated by the numerical results. The numerical relevance 
of this violation (i.e. possible systematic errors) as well as the its physical 
implications are an interesting open problem. 

The contribution by E. Meggiolaro covers in more detail these points. 

4. 3 Increasing density 

Till two or three years ago, we thought that asymptotic freedom was the main 
physical agent behind the pattern of chiral symmetries at high density. Now it 
has been recognized that, at least at zero temperature, the main mechanisms 
are instability at the Fermi surface leading to color superconductivity The 
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rc'vic;w by NarduUi and Sannino cover these points which, unfortunately, are 
not yet amenable to a lattice study. At non-zero density the lattice approach 
(see contribution by M. D'Elia for developments presented at this meeting) 
is still limited to rather high temperature, where color superconductivity is 
most likely lost anyway. 

5 QCD deconfinement transition (rriq = oo) 

When vriq = oo quarks are static and do not contribute to the dynamics: 
hence, the dynamic of the system is driven by gluons alone, i.e. we are dealing 
with a purely Yang-Mills model: 

S = F^^F^"' (22) 

In addition to the local gauge symmetry, the action enjoys the global symme- 
try associated with the center of the group, Z{Nc). The order parameter is 
the Polyakov loop P 

In practice, P is the cost of a static source violating the Z{Nc) global sym- 
metry. 

The interquark potential V(R,T) (R is the distance, T is the temperature) 

is 

^-V(R,T)/T p(0)pt(i|) > (24) 

Confinement can then be read off the behavior of the interquark potential at 
large distance. When V{R) cx <jR it would cost an infinite amount of energy 
to pull two quarks infinitely apart. Above a certain critical temperature V{R) 
becomes constant at large distance: i.e. the string tension is zero, confinement 
is lost. The implication of this is that |Pp — V^(cxd, T) is zero in the confining 
phase, different from zero otherwise. P plays then a double role, being the 
order parameter of the center symmetry, and an indicator of confinement. 
We learn that in Yang Mills models there is a natural connection between 
confinement and realization of the Z{Nc) symmetry. Hence, the confinement 
/ deconfinement transition in Yang Mills systems is amenable to a symmetry 
description. By applying now the same dimensional reduction argument as 
above, we conclude that the Universality class expected of the three color 
model is the same as the one of a three dimensional model with Z{3) global 
symmetry: Indeed, work by the Columbia, Tsukuba and APE group in the 
mid 80's - see again Nicola Cabibbo's talk for comments on this - found that 
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the transition turns out to be 'almost' second order, i.e. very weakly first 
order, like the 3d three state Potts model. 

The same reasoning tells us that the two color model is in the universality 
class of the three dimensional Z{2) (Ising) model. This predictiaa has been 
checked with a remarkable precision by Engels and coUaboratorsllS, and it is 
a spectacular confirmation of the general idea of universality and dimensional 
reduction. 

6 Summary and Open Questions for the QCD High T 
Transition 

What do we know about the real world: two (nearly) massless quark mq < < 
Kqcd, and one more heavy? 

We can approach then the 'real' world from two sides, either decreasing 
the mass from infinity, or increasing the quark mass from zero. 

6. 1 Approaching the physical point from infinite mass 

Remember that in the infinite mass limit QCD reduces to the pure gauge 
(Yang Mills) model. Yang Mills systems have a deconfining transition associ- 
ated with the realization of the global Z{Nc) symmetry. This places the sys- 
tem in the Ising 3d universality class for two colors, and makes the transition 
weakly first order (near second, in fact) for three colors. General universality 
arguments are perfectly fulfilled by the deconfining transition. 

The Z{Nc) symmetry is broken by the kinetic term of the action when 
the quarks are dynamic (m^ < oo) : this particular symmetry description of 
deconfinement only holds for infinite quark mass. When light quarks enter the 
game, the global Z{Nc) symmetry observed at infinite mass is lost, and the 
simple description of confinement in terms of such symmetry is not possible 
any more. It should however remain true that color forces at large distance 
should decrease with temperature: the main mechanism, already at work at 
T = 0, is the recombination of an (heavy) quark and antiquark with pairs 
generated by the vacuum: QQ — > qQ + qQ. At high temperature it becomes 
easier to produce light qq pairs from the vacuum, hence it is easier to 'break' 
the color string between an (heavy) quark and antiquark QQ. In other words, 
we expect enhanced screening of the color forces, which should be sharp at 
a phase transition (or rapid crossover). It is however worth mentioning that, 
even if the string 'breaks' bound states might well survive giving rise to a 
complicated, non-perturbative dynamics above the critical temperature. The 
physical scale of these phenomena is the larger physical scale in the system, 
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i.e. the pion radius. 

6.2 Approaching the physical point from zero mass 

For zero bare mass the phase transition is chiral. For three colors, two flavors 
is second order with Tc ~ 170 MeV. The prediction from dimensional reduc- 
tion + universality -0(4) exponenents- is compatible with the data, but the 
agreement is not perfect. 

If the agreement were confirmed, that would be an argument in favor 
of the non-restoration of the 1/^(1) symmetry at the transition, which is also 
suggested by the behavior of the masses spectrum. Remember in fact that the 
chiral partner of the pion is the /o, which is in turn degenerate with the scalar 
oq with Ua{^) is realized. All in all, C^a(1) non-restoration across the chiral 
transition corresponds to ~ m^g ^ ruag which is the pattern observed in 
lattice calculations (once more I refer to meggiolaro's talk). 

The transition with three (massless) fiavor turns out to be first order. 
The question is than as to whether the strange quark should be considered 
'light' or 'heavy'. In general, the real world will be somewhere in between 
two and three light flavor, and to really investigate the nature of the physical 
phase transition in QCD one should work as close as possible to the realistic 
value of the quark masses, which is a very demanding numerical task. 

6.3 What do we know on the real QCD phase diagram 

Among the most prominent open questions, there is of course the behavior of 
'real' QCD, with two light flavor, and a third one of the order of Aqcd, so 
how and when exactly the Nf = 2 scenario morfs with the Nf — 3? Also, 
why is much smaller that the pure gauge deconfining transition? 

At a theoretical level the question is if it is possible to give an unified 
description of the two transitions, chiral and deconfining. This question is 
currently under active investigation: recent work suggests that a symmetry 
analysis of the deconfining transition can be extended also to theories with dy- 
namical fermions. The physical argument is rooted in a duality transformation 
which allows the identification of magnetic monopoles as agent of deconfine- 
ment. The xirder parameter for deconfimenent would that be the monopole 
condensate El. An alternative approach uses percolation as the common agent 
driving chiral and confining transitions □. 

One unifying description and perhaps the most dramatic evidence of a 
phase transition away from the 'simple' limits ruq — and — oo comes from 
the equation of state: for T ~ 180 Mev we observe, from lattice calculations, 
a sharp increase of the internal energy: the behavior of the internal energy is 
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a direct probe o£j.he number of degrees of freedom, and indicates quark and 
gluon liberation 113. _ 

Finally, the work 113 arrives at an interesting picture of the phase diagram 
of QCD by combining symmetry analysis and phenomenological considera- 
tion. Particularly interesting is the prediction of an endpoint of a first order 
line stemming from zero temperature chiral transition at finite density, which 
should be experimentally observable. 



7 Methods for QCD at finite temperature and density 

Here we will concern ourselves with computational schemes for QCD. The 
methods described here are those used to obtain the results reviewed above 
on the chiral transition, deconfinement and equation of state in QCD. We 
will give some details on these methods, we will explain why they are not 
immediately applicable at finite density and we will close up with a brief 
assessment of the current situation for finite density QCD. 

The question is how to estimate the physical observables < O >= 
TrOp/Z starting from the representation of the partition function 



Z{n,T)^j dty e-^«+'^(^+"+''^°)'^dV^dV'dC/ (25) 

The need for two integrations (over bosons and over fermions) lends itself 
naturally to two different paths: either integrate gluons first, or fermions. 

In the first case (if gluons are integrated out first) we end up with a purely 
effective fermionic model: 

z(T,Ai,^,^,c/)^z(r,/i,v;,^) (26) 

As the integration over gluons cannot be done exactly |^, such effective models 
are often build with the help of a symmetries' analysis for QCD. The gauge 
fields enters the the game under the guise of coefficients for such models. 
Very important examples of this approach include instanton models, chiral 
perturbation theory, four fermion models such as those discussed by Nardulli 
and Sannino at this meeting. 

The other way takes advantage of the bilinear nature of S in the fermionic 
fields, yielding the exact expression: 

Z{T,fi,U)^ f dUe^^^^~^°s{detM)) (27) 



There indeed a a systematic approach to QCD based on the strong coupHng expansion 
which makes this integrations easy, but this is not in the scope of this introduction 
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The above integral needs being regulated, and scheme for calculating it has 
to be devised. Both tasks are accomplished within the lattice approach. 

7.1 Lattice field theory at T,ii ^ 

Temperature comes for free on a lattice: the lattice has a finite extent Nta, 
hence temperature is given hy T ~ 1/Nta. The discretization can be carried 
out in complete analogy to T = 0, and most of the techniques developed there 
(see e.g. the recent review for an introduction to lattice field theory at zero 
temperature) apply at finite temperature as well. 

A finite density of baryons /is Jq |^ poses instead specific problems. 

How to discretize fiipjo^ 

Recall first the 'natural' discretization of the matter fields and their 
derivatives d^4>(x) on a regular lattice with spacing a: 

(t>LATT{ni,n2, na, n^) = 4>{nia, n2a, n^a, n^a) 
Af^(l)LATT{ni,n2, na, n^) = (0(nia, (n^ + l)a, nsa, n^a) - 

0(nia, (n^a, n^a, n4a))/2a 

This, for instance, is the correct prescription for the chiral condensate: 
rmpip — > rnipLATT'4'LATT but not, as we will see in a moment, for the baryon 
density! fiipjoip P^'4'lattIo'^latt is not the correct lattice form , i.e. The 
naive discretization is not adequate for baryon density. 

Let us then consider free fcrmions in the continuum 

_ _ _ 

Jo 

The internal energy e = ^^InZ ^ ^ d^pj^JS0^^^^_^^ after sub- 

4 

racting the vacuum energy, is finite at T = 0: liniT^o^ = a-nd gives the 
expect result 

By use of the naive discretization, e would instead diverge in the contin- 
uum (a ^ 0): i = i'xjfiipx+fj.a + rwipxipx + fJ''4'xJo'4'xi oc ^ ^a^o 00 

The solution 0: /i is an external field 

Note the analogy: ?/)7pA^'0 < > ifi'ip'^Qip. It shows us that /i looks like an 

external field in the time direction. But we know how to 'put' external field 

Remember again that = Lq + M'^Oi -^0 = V'7oV'i N — N = J Jq 
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on a lattice: they live on the lattice links: for instance, in electrodynamics 
A ^ 9 = e^^ (and, for the gauge fields, this implements in a natural and 
elegant way gauge invariance) 

These considerations suggest how to put finite density on a lattice: 

= V^.7oe''"^,+6 - V^x+o7oe-'"'Vx (28) 

Indeed, this turns out to be the correct prescription: unphysical divergences 
disappear and the continuum limit is reproduced. 

There is also a very expressive physical interpretation: as we can see, 
forward propagation is encouraged and backward propagation is suppressed: 
we are indeed inducing an asymmetry particles-antiparticles. 

The 0-th component of the current jo counts indeed the differences be- 
tween backward and forward propagation: 

Jo = -d^,L = V^xloe'^Vc.+o + V^x+67oe-''"Vx (29) 
and reproduces the correct continuum limit. 

Chemical potential and boundary conditions 

By use of an unitary transformation of the fields it is possible to re-express 
L as 



L(/x) = m (30) 

with boundary conditions 

(j){x + NT) =e^''^''(j){x) (31) 

ip{x + Nt) = -e'^^'^tpix) (32) 

It is of some interest to consider the effect of the chemical potential on 
the baryonic propagators in the two cases - when the chemical potential is 
included into the Action, and when instead affects only the boundary condi- 
tions. In the first case, (at zero temperature, and ignoring feedbacks) a term 
exp{—NcH) multiplies the baryonic propagator: this produces an apparent 
decrease of the baryon mass 

niB = niB — 3/z (33) 

i.e. the baryon becomes massless at n = = ttib/S. In the same situation, 
the chiral condensate remains constant till fXc, and then suddenly drops to 
zero. The behavior of the chiral condensate would be suggestive of a strong 
first order transition, while the behavior of the baryon mass would suggest a 
second order transition with u = 1\ 
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In the; other fornnilation (when the chemical potential only affects the 
boundary conditions), at zero temperature the baryon mass is constant till /Uc 
: this is consistent with the behavior of the chiral condensate, and with the 
physical intuition that nothing should happen till the Fermi level is reached. 
Of course, nothing is wrong and the apparent contradiction is resolved by 
noticing that the apparent decrease of the baryon mass below /Uc merely re- 
flects the change of the reference energy. The relationships above provide a 
link between the two pictures. 

QCD at nonzero T and /i at a glance 



The continuum formulation: 

L = Lym + ipiijixD^ + m)il) 
+ iJ.ip'yo'tp 

(34) 

/U is explicitly included via the coupling to /i Jq , and the temperature is the 
reciprocal of the imaginary time. 
On a lattice: 

3 



^ appears as a link term, and the temperature is again the reciprocal of 
the imaginary time. The two formulations coincide in the limit a — > 0. 

7.2 From the formulation to the results 
Let us write again: 

S = Sym{U) + i)M{U)tlj (35) 
By taking advantage of S bilinearity in the fields V, V' we can write 

^ = J e-'^^'^^^UetM{U)dU (36) 

It is convenient to define an 'effective' Action 

Seff{U) = Sym{U) - ln{detM{U)) (37) 
Averages of purely gluonic observables can be expressed as 

< fiU) >= Z-^ j dUe-^'ff^^^f{U) (38) 
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while fermion bilinears can be evaluated with the help of Grassman algebra 
< qq >= Z-^ I dUe-^-"^^'>DetM{Uy^ (39) 



In conclusion, if we know how to treat / e~^^f^'^^^dU we have access to all 
of the gluonic observables, chiral condensate, meson and baryons propagators 
(masses, decay constant, etc.) etc. much in the same way as at T = 0. 

In practical numerical works lattice discretization is combined with im- 
portance sampling: a configuration of gauge fields [U] is a 'point' in a mul- 
tidimensional integration space. A Markov chain of points is then created 
according to the prescription: P([J7]) oc e^^'^f'^^^^^ Expectation values are 
then given by simple averages: < O >= liniff^oo^/N^-^^ 0{U). 

The prescription above relies on importance sampling: Seff{\U]) must 
be positive ^ In QCD (eq. above) M'^{p) — — Af(— /i) : importance sam- 
pling, hence MonteCarlo evaluation of physical observables, works for purely 
imaginary chemical potential, including, of course, /z = 0. 

For QCD with two color, however, as well as for several fermionic models, 
the positivity condition is met by the Action at non zero /i. All in all, we have 
a rather elaborate pattern of warnings and possibilities. 

7.3 Results and possibilities for QCD at fi ^ 

I would classify the many questions one can ask along two main lines: 

Firstly, what can be done in order to learn about thegeneral behavior at 

nonzero baryon density? 

Secondly, what can be done in order to learn about real QCD at nonzero 

baryon density? 

I would say that to learn about the effects of a finite density of baryons, we 
can study two color QCD : here, not only one can study fermionic observables, 
but it is also possible to observe the effects of a density of baryons on the gauge 
fields. This is indeed an active field of research, and I refer to the proceedings 
of the Lattice meetings for details. 

How about real QCD? there is a growing consensus that at least the high 
temperature, small density region is accessible: by increasing temperature it 
is easier to fluctuate light baryons, hence to explore a region with a moder- 
ate baryon density. At least three different, complementary approaches stem 
from this consideration:direct evaluation of derivatives at /z = 0, reweighting, 

'^If this is not the case, one might think of using instead: this has to be used with care 
because it introduces extra degrees of freedom which might well be dangerous. 
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analytic continuation from purely imaginary chemical potential M. This last 
approach has been discussed by M. D'Elia at this meeting. 

I hope to have managed at least to describe some of the aspects of a 
rapidly evolving field, and to convey a feeling of its many reasons of interest. A 
discussion of the recent results is not in the scope of this cursory introducticp: 
for this, and for a complete set of references, I refer to the latest reviews Ej 
(not yet available at the time of the Frascati meeting). 
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